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ABSTRACT 


A  location  study  is  made  of  28  large  underground 
explosions  detonated  in  the  northern  area  of  the 'Nevada 
Test  Site  (NTS).  Recording  networks  were  comprised  of 
between  9  and  49  teleseismic  stations  having  two-  to  three- 
quadrant  distributions. 

Errors  of  locations  obtained  without  applying  travel¬ 
time  anomalies  (relative  residuals),  and  with  depths 
restrained  to  the  known  values,  average  about  7  km  but  are 
as  large  as  20  km.  With  anomalies,  the  errors  are  consist¬ 
ently  2.  5-3.0  km.  The  size  of  the  area  at  NTS  across  which 
the  anomalies  are  valid  is  at  least  70  km  by  25  km. 

Depth  errors  average  70  km  without  anomalies  and  15  km 
with  anomalies.  Without  the  anomalies,  a  linear  relation¬ 
ship  is  observed  between  the  least-squares  standard 
deviation,  a, of  time  errors  of  the  solution  obtained  when 
the  depth  is  restrained  to  its  true  value  and  the  depth 
errors,  dz,  of  the  corresponding  unrestrained  solution: 
dz(km)=75a  (sec). 

By  deliberately  mislocating  a  calibration  event 
approximately  140  km,  it  is  shown  that  relative  accuracy 
(precision)  remains  at  about  2, 5-3.0  km. 

If  a  constant  network  is  used  to  locate  a  set  of 
events,  the  "bias"  is  a  consistent  and  unique  function 
of  that  network.  The  lack  of  a  common  bias  among  networks 
with  similar  azimuthal  distributions  definitely  eliminates 
the  source  region  as  the  principal  cause  of  the  anomalies. 
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Alternatively,  the  anomalies  may  be  attributed  to  slight 
lateral  and  vertical  inhomogeneities  within  the  mantle 
between  the  source  and  receiver,  the  effects  of  which  are 
integrated  along  the  entire  path. 

Furthermore,  when  constant  networks  are  used,  the 
relative  locations  obtained  without  anomalies  are  identical 
to  those  obtained  with  anomalies,  except  for  a  bias  trans¬ 
lation  appropriate  for  that  network. 

Finally,  it  is  shown  that  station  anomalies  determined 
from  explosions  occurring  in  regions  other  than  NTS  do  not 
agree  with  those  at  NTS.  It  is  demonstrated  that,  in 

general,  it  is  better  to  apply  no  anomalies  at  all  rather 
than  the  wrong  ones. 
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INTRODUCTION 


In  earlier  papers  (Chiburis,  1968b,  Evernden,  1969), 
it  was  demonstrated  that  a  significant  improvement  in 
accuracy  can  be  achieved  when  locating  nuclear  explosions 
with  networks  of  limited  azimuth  aperture  (single  quadrant) 
and  limited  numbers  of  stations,  if  pre-determined  travel¬ 
time  anomalies1  are  used.  The  accuracy  reported  for  some 
seventeen  well  recorded  explosions  in  the  Nevada  Test  Site 
area  (NTS)  was  2. 5-3.0  km  when  the  solutions  were  restrained 
to  the  surface.  The  objectives  of  the  present  study  are 
to  investigate  in  more  detail  the  stability  and  applica¬ 
bility  of  the  anomalies  across  and  beyond  NTS,  the  effects 
of  location  bias  on  anomaly  solutions,  unrestrained  depth 
solutions,  and  constant  and  variable  network  locations. 

It  is  important  to  understand  the  definition  of  the 
term  "bias"  In  terms  of  the  particular  network  used  for 
any  location,  and  to  understand  that  the  results  of  any 
solution  are  relatively  unpredictable  when  no  heed  is  taken 
of  the  anomalies.  In  this  report,  "location  bias"  is 
defined  as  that  error  obtained  when  anomalous  travel  times 
are  used  for  locating  teleseismic  explosions,  regardless  of 
the  causes  of  the  anomalies.  The  travel  times  are  termed 


An  anomaly  is  defined  as  the  usual  station  residual  but 
relative  to  another  station  in  the  network;  that  is,  if 
the  residual  at  station  i  is  R.  =  observed  time  -  computed 
time,  then  the  anomaly  at  station  i  relative  to  station  j 


anomalous  in  that  they  do  not  conform  to  those  of  an  average 
earth  model.  It  is  generally  agreed  that  the  anomalies  are 
caused  by  unknown  inhomogeneous  effects  of  source  region 
geology,  of  the  crust  and  upper  mantle  in  the  vicinity  of 
the  recording  station,  and  of  the  total  travel  path  in  the 
deeper  mantle.  Disagreement  arises  as  to  the  proportionate 
contribution  of  each  of  these  effects.  At  the  present  time, 
no  technique  is  known  which  can  unambiguously  isolate  the 
three  possible  causes  or  quantify  their  respective  contri¬ 
butions.  Although  the  exact  causes  of  the  anomalies  are 
not  considered,  it  will  be  shown  that  their  effects  on 
£el.^.t:‘-ve  location  accuracy  using  a  common  network  are  negli¬ 
gible  (within  the  limits  of  seisinograni  timing  error),  and 
that  if  one  does  not  compensate  for  the  existence  and  vari¬ 
ability  of  anomalies  for  different  stations  and  regions, 
the  calculated  locations  may  be  highly  variable. 
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DESCRIPTION  OF  THE  DATA 


A  total  of  28  large  (M=4.8-6.4)  explosions  detonated 
within  the  northern  portion  of  NTS  (Figure  1)  was  selected 
for  analysis  on  the  basis  of  the  number  of  recording  stations 
and  the  quality  of  the  signals  received.  Information  pertinent 
to  these  events  is  given  in  Table  I,  which  also  includes  the 
number  and  type  of  recording  stations  used  for  location 
purposes.  All  of  the  selected  stations,  a  total  of  65,  were 
teleseismic  (>  1600  km)  to  the  explosions.  Table  II  lists 
the  stations,  their  type,  distances  and  azimuths  from  the 
event  BILBY,  and  whether  they  recorded  a  particular  event. 

All  time  readings  were  made  by  analysts  at  the  Seismic  Data 
Laboratory  ard  are  believed  accurate  to  within  0.1  sec  for 
most  of  the  LRSM  stations,  within  0.5  sec  for  the  USC&GS 
stations,  and  better  than  0.1  sec  for  the  one  VELA  Observ¬ 
atory  (CPSO).  The  travel-time  tables  used  for  the  locations 
are  those  of  Herrin  (1968),  unless  clearly  specified  other¬ 
wise. 
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TABLE  II 

NETWORK/EVENT  COMBINATIONS 


LOCATION  PROCEDURES  AND  RESULTS  USING  NO  ANOMALIES 


To  determine  the  improvement  in  accuracy  it  is  first 
necessary  to  locate  all  explosions  with  every  available 
recording  station,  to  restrain  the  depths  to  the  known 
values,  and  to  use  no  anomaly  corrections.  The 
number  of  stations  recording  each  of  the  events  varies 
from  9  to  49;  the  azimuth  aperture,  or  quadrantal  coverage, 
varies  from  123°  to  245°  (two  to  three  quadrant  distri¬ 
bution).  The  location  results  are  shown  in  Figure  2  both 
in  plan  view  and  normalized  to  a  common  origin.  The 
numerals  adjacent  to  the  vectors  in  Figure  2  refer  to  the 
event  identification  number  as  given  in  Figure  1.  In  order 
to  note  any  possible  variations  in  location  patterns,  the 
southern  area  of  NTS  was  split  into  two  subregions,  I  and 
TT*  The  location  errors  for  the  28  events  vary  from 
0.60  km  to  18.89  km  (with  an  average  of  6.65  km)  and  the 
directions  of  the  shifts  are  nearly  random,  although  there 
is  some  suggestion  that  in  subregion  I  the  events  shift 
northerly  and  in  subregion  II  southwesterly.  Since  the 
networks  are  well-distributed,  and  the  s ignal-to-noise 
values  are  high,  these  results  are  indicative  of  the  smallest 
errors  one  can  expect  to  achieve  for  the  NTS  region  without 
empirical  corrections. 

At  this  point,  it  is  of  interest  to  compare  the 
preceding  results,  obtained  with  the  Herrin  1968  tables, 

with  those  obtained  with  the  Jeffreys-Bullen  (J-B)  tablesT - 

The  d— B  results  are  shown  in  Figure  3  in  which  the  errors 
vary  from  2.69  km  to  21.11  km  (with  an  average  of  7.66  km). 
The  difference  in  the  average  error  between  the  two  tables 
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io  negligible,  although  individual  events  have  quite 
different  shift  magnitudes,  and  the  normalized-origin  plot 
is  considerably  changed  from  that  in  Figure  2:  the  events 
in  subregion  I  shift  nearly  in  one  quadrant  (north-north- 
easterly)  and  those  in  subregion  II  the  same  except  for 
one  event.  The  location  results  for  both  travel-time 
tables  are  given  in  Table  III,  which  also  includes  the 
standard  deviations  of  errors  obtained  from  the  least- 
squares  solutions.  The  values  of  the  standard  deviations, 
which  are  general  indicators  of  the  goodness-of-fit  to 
the  least-squares  regression,  are  of  the  same  size  as 
those  obtained  when  locating  events  anywhere  with  any 
network  (e.g.,  those  reported  by  the  USCSGS  on  PDE  cards). 
Therefore,  the  times  read  in  this  study  are  not  unusually 
good  regarding  the  least-squares  fits. 

Generally,  for  the  NTS  region,  the  results  indicate 
that  when  anomalies  are  not  used,  1)  restrained  location 
errors  are  moderate  on  the  average  (about  7  km)  but  can 
be  large  (about  20  km);  2)  no  consistent  location  bias 
exists  which  can  be  reasonably  attributed  to  NTS  source 
geology  alone;  and  3)  although  both  the  Herrin  and  J-B 
tables  yield  about  the  same  magnitude  of  average  error, 
the  patterns  of  locations  are  a  function  of  the  particular 
travel-time  table  used. 
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LOCATION  PROCEDURES  AND  RESULTS  USING  CALIBRATION  ANOMALIES 


Following  the  procedures  described  by  Chiburis  (1968b), 
relative  travel-time  anomalies  were  determined  from  the 
events  BILBY,  TAN,  COMMODORE,  and  CORDUROY.  Those  stations 
recording  more  than  one  of  these  "calibration  events"  were 
accordingly  assigned  an  anomaly  based  on  the  average  of 
from  two  to  all  four  events.  The  resulting  anomalies, 
relative  to  station  RK-ON  are  presented  in  Table  IV. 

Assuming  that  the  anomalies  from  these  four  events  effec¬ 
tively  correct  the  inadequately-known  -earth  model  appropriate 
to  NTS,  the  remaining  24  events vcan  be  located  using,  as 
before,  all  available  stations  for  each  event.  The  location 
shifts  are  plotted  in  Figure  4  and  the  results  given  in 
Table  V.  Using  anomalies,  the  location  errors  now  range 
from  0.52  km  to  6.75  km  with  an  average  of  2.81  km,  an 
improvement  factor  of  about  three  over  the  results  obtained 
without  anomalies.  The  corresponding  standard  deviations 
are  significantly ’ reduced  (a  necessary  but  not  sufficient 
condition  for  the  anomalies  to  be  appropriate  for  the 
NTS  region) . 

Although  for  some  purposes  an  improvement  factor  of 
three  might  not  appear  to  warrant  the  general  use  of  anomalies, 
there  are  several  facts  to  consider:  1)  the  locations 
without  anomalies  were  made  with  an  average  of  26  recording 
stations  which  were  well-distributed  azimuthally  such  that 
the  networks  were  quite  stable  (in  fact,  even  nearly  ideal 
network  distributions  do  not  guarantee  more  accurate 
locations,  as  Lambert  et  al  (1969)  report  a  20  km  error 
for  LONG  SHOT  when  using- 329  stations  at  virtually  all 
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TABLE  IV 


f 


CALIBRATION 

ANOMALIES; 

HERRIN 

1 96B  TRAVEL 

-TIME  TABLES 

Station 

(VELA) 

Anoma ly 
-(Sec) 

Source 

Code* 

Station 
( C4GS ) 

Anoma ly 
(Sec) 

Source 

Code* 

Station 

(C&GS) 

Anoma  ly 
(Sec) 

AD-  IS 

0.16 

D 

AAM 

0.03 

C 

MAL 

0.  58 

AX2AL 

1.68 

D 

AKU 

1.04 

D 

MAT 

-0.21 

BE-FL 

1.24 

C 

ANT 

1 .  16 

C 

NOR 

-0. 4B 

BL-WV 

0.56 

D 

ARE 

1  .36 

B 

NNA 

0.65 

BR-PA 

0.71 

D 

ATL 

0.86 

C 

NUR 

-0.31 

CPSO 

0.86 

B 

BEC 

e 

0.35 

A 

OGD 

0.22 

OH-NY 

r°.02 

D 

BLA 

0.95 

C 

OXF 

1  .  73 

EB-MT 

0.09 

0 

CAR 

0.63 

A 

PEL 

0.12 

EN-MO 

0.68 

0 

COL 

0.91 

B 

PTO 

-0.54 

EU-AL 

1.98 

D 

COP 

0 .  l  ;■ 

D 

SCP 

0.19 

EU2AL 

1.98 

•  ** 

CMC 

-U .  G 1 

B 

SEO 

0.01 

GG-GR 

0.96 

D 

ESK 

-0.77 

D 

SHA 

3.13 

HN-ME 

0.56 

A 

FLO 

1.45 

C 

SHK 

-0.45 

LV-LA 

1.95 

D 

GEO 

0.29 

B 

SJG 

0.72 

NP-NT 

M.0| 

B 

GOH 

0.60 

B 

STU 

0.46 

00-NW 

-0.04 

0 

GIE 

1.70 

D- 

TOL 

0.26 

PG-BC 

2.45 

0 

GUA 

-1.44 

C 

TRN 

-2.16 

PZ-PR 

1.22 

D 

KEV 

-1.19 

C 

VAL 

-0.66 

RK-ON 

0 

** 

KIP 

1.64 

B 

WES 

0.60 

SI-BC 

2.90 

C 

KON 

0.26 

D 

SV2QB 

-0.42 

*** 

KTG 

0.05 

B 

SV3QB 

-0.42 

B 

LPB 

0.12 

B 

WH2YK 

0.75 

D 

LPS 

1.63 

B 

Source 

Code* 

C 

B 

B 

B 

B 

B 

A 

B 

B 

A 

C 

C 

D 

B 

A 

A 

A 

D 

D 


** 

*** 


Code  =  A  Anomalies  are  an  average  of  BILBY,  TAN,  CORDUROY,  COMMODORE. 

Code  =  B  Anomalies  are  an  average  of  three  of  the  above. 

Code  =  C  Anomalies  are  an  average  of  two  of  the  above. 

Code  ■  D  Anomalies  were  determined  from  only  one  of  the  above. 

All  anomalies  are  referenced  to  station  RK-ON. 

Used  anomaly  determined  from  EU-AL  for  EU2AL  and  SV2QB  for  SV3QB. 
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TABLE  V 


LOCATION  RESULTS  USING  CALIBRATION  ANOMALIES 


Standard 


Event 

Number 
Sta  ti ons 

Azi muth 
Aperture 

Shift 

(km) 

Di rection 
of  Shift 
(Deg) 

Devi ati on 
of  Errors 
(Sec) 

GREELEY 

49 

194 

2.86 

51 

0.40 

CHARTREUSE 

28 

194 

3.60 

225 

0.31 

SCOTCH 

26 

178 

4.15 

193 

0.38  * 

KNICKERBOCKER 

22 

195 

1.94 

196 

0.45 

DURYEA 

21 

178 

1.83 

149 

0.35 

BOXCAR 

17 

232 

3.  39 

116 

0.30 

REX 

16 

157 

2.57 

197 

0.40 

BENHAM 

11 

201 

0.94 

192 

0.32 

DUMONT 

40 

194 

1.84 

275 

0.22 

PILEDRIVER 

38 

245 

2.03 

179 

0.33 

PIRANHA 

36 

234 

1.37 

210 

0.41 

BRONZE 

35 

238 

0.52 

136 

0.36 

AGILE 

35 

238 

2.36 

213 

0.45 

WAGTAIL 

29 

188 

2.30 

250 

0.30 

FORE 

26 

138 

5.85 

248 

0.43 

BUFF 

25 

238 

1.85 

91 

0.35 

NASH 

22 

185 

4.17 

43 

0.39 

TURF 

21 

157 

1.06 

41 

0.38 

CUP 

20 

234 

3.45 

222 

0.29 

KLICKITAT 

18 

157 

2.99 

236 

0.36 

BOURBON 

16 

185 

1.93 

183 

0.43 

AUK 

11 

155 

6.75 

44 

0.35 

PINSTRIPE 

11 

123 

6.14 

80 

0.25 

PAR 

9 

182 

1.55 

47 

0.30 

AVERAGE  24  2.81 


azimuths  and  distances);  2)  in  a  previous  paper  (Chiburis, 
1968b),  where  the  networks  were  approximately  single-quadrant 
and  had  fewer  recording  stations,  a  more  common  situation, 
the  location  results  without  anomalies  were  20-25  km  in  error 
and  3)  in  both  the  previous  and  present  studies,  the  location 
accuracies  obtained  using  anomalies  are  essentially  the  same 
(2. 5-3.0  km)  regardless  of  network  stability.  Therefore,  all 
locations  fall  in  a  small  area  around  the  true  location  with 
the  use  of  anomalies,  whereas  otherwise  the  solutions  cover 
a  large  area,  are  mislocated  (in  an  absolute  sense  if  the 
true  locations  happen  to  be  known,  as  for  explosions),  and 
are  prone  to  misinterpretation  regarding  their  relationship 
to  assumed  earth  models,  to  tectonic  or  explosion  patterns, 
and  to  source  region  geophysics. 


USE  OF  AVERAGE  ANOMALIES 

Seismograms  of  any  event  have  some  random  reading 
errors  associated  with  them  which  can  affect  the  location 
of  subsequent  events  when  using  the  measured  anomalies  as 
corrections.  The  location  results  in  Table  V  obtained 
with  anomalies  indicate  the  expected  accuracy  when  one  to 
four  events  are  available  for  calibration  purposes.  To 
minimize  these  random  errors,  an  average  anomaly  can  be 
determined  for  each  station  from  as  many  events  within 
the  same  anomaly  region  as  possible.  This  process  will 
give  the  best  correction  for  that  region.  Table  VI  gives 
such  corrections  and  Figure  5  shows  the  location  shifts 
of  the  28  events  obtained  by  using  the  average  anomalies. 
The  results  are  tabulated  in  Table  VII,  where  the  average 
error  is  seen  to  be  2.14  km,  which  represents  perhaps  the 
optimum  accuracy  for  teleseismic  location  at  NTS  when  the 
standard  deviation  of  reading  errors  is  on  the  order  of 
0.25  sec . 


TABLE  VI 


AVERAGE  ANOMALIES  FOR  NTS;  HERRIN  1968  TRAVEL-TIME  TABLES 


Station 

(VELA) 

Number 

Anoma ly 
(Sec) 

of 

Even  ts  * 

S  tation 
(C4GS  ) 

AD- 1 S 

0.01 

6 

AAM 

AX2AL 

1.85 

8 

AKU 

BE-FL 

1.35 

6 

ANT 

BL-WV 

0.75 

7 

ARE 

br-pa 

1.00 

8 

ATL 

CPSO 

0.93 

22 

BEC 

DH-NY 

0.17 

7 

BLA 

EB-MT 

0.26 

4 

CAR 

EN-MO 

1.06 

3 

COL 

EU-AL 

2.01 

2 

COP 

EU2AL 

2.01 

*** 

CMC 

GG-GR 

0.69 

4 

ESK 

HN-ME 

0.64 

25 

FLO 

lv-la 

1.96 

2 

GEO 

NP-NT 

1.05 

24 

GDH 

OO-NW 

-0.07 

4 

GTE 

PG-BC 

3.06 

11 

GUA 

PZ-PR 

1.56 

2 

KEV 

RK-ON 

0 

** 

KIP 

SI-BC 

3.06 

1C 

KON 

SV2QB 

-0.45 

4 

KTG 

SV3QB 

-0.47 

15 

LPB 

WH2YK 

1.25 

3 

LPS 

Anomaly 

_(.s.ecj_ 

Number 

of 

Even  ts  * 

Stati on 
( C&GS ) 

Anomaly 

i.Sec) 

Number 

of 

Events* 

0.51 

16 

MAL 

0.55 

1 1 

1.08 

7 

MAT 

-0.18 

10 

1.12 

1 1 

NOR 

-0.70 

15 

1.30 

20 

NNA 

0.37 

7 

1.00 

19 

NUR 

-0.28 

14 

0.65 

9 

OGD 

0.50 

1 1 

1.38 

13 

OXF 

1.99 

24 

0.83 

19 

PEL 

-0.12 

8 

1.34 

26 

PTO 

-0.74 

8 

0.13 

2 

SCP 

0.43 

16 

0.24' 

18 

SEO 

0.42 

13 

-0.77 

7 

SHA 

3.36 

8 

1.63 

20 

SHK 

-0.10 

5 

0.45 

10 

SJG 

0.71 

18 

0.72 

11 

STU 

0.46 

13 

2.00 

3 

TOL 

0.17 

15 

-1.43 

6 

TRN 

-2.04 

13 

-1.13 

7 

VAL 

-0.38 

3 

1.89 

11 

WES 

0.54 

2 

0.12 

13 

0.17 

16 

***  ,,  an°malles  are  referenced  to  station  RK-ON. 
sed  the  anomaly  determined  for  EU-AL. 
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TABLE  VII 

LOCATION  RESULTS  USING  AVERAGE  ANOMALIES 


Event 

Number 

Stations 

Azimuth 

Aperture 

(Deg) 

Shift 

(km) 

Di rection 
of  Shift 
(Deg) 

Standard 
Deviation 
of  Errors 
(Sec) 

GREELEY 

49 

194 

3.66 

37 

0.28 

CHARTREUSE 

28 

194 

3.45 

235 

0.26 

SCOTCH 

26 

178 

1.89 

199 

0.32 

KNICKERBOCKER 

22 

195 

0.61 

276 

0.32 

DURYEA 

21 

178 

0.95 

200 

0.28 

BOXCAR 

17 

232 

1.91 

92 

0.17 

REX 

16 

157 

0.45 

225 

0.32 

BENHAM 

11 

201 

0.72 

280 

0.24 

DUMONT 

40 

194 

2.00 

275 

0.17 

PILEDRIVER 

38 

245 

1.19 

177 

0.34 

PIRANHA 

36 

234 

1.22 

237 

0.34 

BRONZE 

35 

238 

0.78 

58 

0.37 

AGILE 

35 

238 

1.89 

251 

0.37 

WAGTAIL 

29 

188 

0.27 

339 

0.24 

FORE 

26 

138 

2.71 

234 

0.34 

BUFF 

25 

238 

2.26 

87 

0.32  1 

NASH 

22 

185 

4.37 

20 

0.34  ( 

TURF 

21 

157 

2.99. 

57 

0.29 

CUP 

20 

234 

3.33  ' 

213 

0.18 

KLICKITAT 

18 

157 

2.17 

230 

0.28 

BOURBON 

16 

185 

1.52 

190 

0.37 

AUK  / 

11 

155 

5.58 

44 

0.27 

PINSTRIPE 

11 

123 

4.47 

72 

0.26 

PAR 

9 

182 

1.84 

55 

0.18 

COMMODORE 

38 

244 

3.13 

23 

0.27 

TAN 

39 

244 

0.78 

242 

0.25 

BILBY 

33 

237 

1.78 

127 

0.20 

CORDUROY 

33 

195 

1.99 

22 

0.25 

AVERAGE 

26 

J 

i 

2.14 

0.28 

/ 

I 
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nts  Anomaly  stability 

I 

Considering  the  location  errors  in  Figures  2  and  4,  there 
appears  to  be  no  significant  difference  between  subregions  I 
and  II  when  using  anomalies  determined  from  events  in  subregion 
I.  The  implication  is/ that  the  region  over  which  the  anomalies 
are  valid  is  at  least  hs  large  as  that  shown  in  Figure  1,  approx¬ 
imately  70  x  25  km.  Although  it  is  known  that  the  anomalies 
vary  from  one  region  to  another  in  an  unpredictable  manner  (e.g., 
Chiburis,  1966a,  b,  and  1968a;  Herrin  and  Taggart,  1968;  Hales 
et  al  (1968);  and  oth/ers),  calibrating  specific  regions  of  interest 
is  possible  if  a  few  i  well-located  shallow  events  are  available. 

This  statement  must  be  qualified  in  that  all  of  the  results  in 
this  study  were  obtained  with  explosions,  and,  hence,  perfectly- 
known  positions,  so  thay  are  undoubtedly  simpler  to  evaluate 
quantitatively  than -those  obtained  with  earthquakes  of  unknown 
position,  especially  depth.  It  is  certain  that  the  anomalies 
change  with  depth  a£  well  as  region,  but  it  is  uncertain  how 
they  change  and  how  best  this  fact  can  be  used  to  advantage. 

Current  studies  are  underway  which  specifically  address  the 
problem  of  travel-time/depth. 
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UNRESTRAINED  DEPTH  SOLUTIONS 

It  is  important  at  this  point  to  investigate  the 
P— wave  solutions  obtained  by  letting  the  depth  parameter 
run  free  in  the  least-squares  scheme.  Using  all  available 
stations,  unrestrained  locations  were  computed  both  with 
and  without  anomalies  determined  from  BILBY,  TAN,  COMMODORE, 
and  CORDUROY.  Table  VIII  gives  the  results  for  the  set  of 
events.  As  indicated  on  the  table,  those  solutions  yielding 
depths  above  the  surface  are  restrained  to  zero  and  the 
epicenter  recomputed.  The  computed  above-surface  values 
shown  in  parentheses  (a  total  of  four)  were  replaced  by 
zeroes  in  forming  the  average.  When  anomalies  are  not  used, 
the  average  depth  error  id  69.4  km  (range:  23.0-121.0  km); 
with  anomalies,  the  average  depth  error  is  reduced  to  14.9 
km  (range.  0.4-42.3  km).  The  unrestrained  epicenter  shifts 
without  anomalies  average  16.75  km  (range:  2.04-58.49  km) 
and  with  anomalies  3.90  km  (range:  0.36-17.36  km).  There¬ 
fore,  for  this  set  of  events,  the  accuracy  is  improved  by 
factors  of  about  five  and  four  in  depth  and  epicenter 
respectively  when  anomalies  are  used. 

These  results  demonstrate  that  if  anomalies  are  determined 
from  teleseismic  explosions,  they  can  be  used  to  teleseismi— 
cally  locate  subsequent  explosions  detonated  within  the  same 
anomaly  region  with  an  epicenter  accuracy  of  about  3  km  and 
a  depth  accuracy  of  about  15  km.  The  same  accuracy  can  be 
expected  for  shallow-focus  (<60  km)  earthquakes  in  the  same 
region,  because  other  studies  (Chiburis  and  Ahner,  1969) 
show  that  the  anomalies  from  explosions  are  approximately 
similar  to  those  from  nearby  shallow  earthquakes. 
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TABLE  VIII 


COMPARISON  OF  DEPTH-TREE  LOCATION  RESULTS  OBTAINED  WITHDUT  ANU  WITH  ANOMALIES 
■  _ Without  Ano  Malles _ _ W  i  th  Anoma  lies 


Actual  Depth 


Event  Name 

Depth 

(km) 

Shift 

(km) 

E  rror 
(  km  ) 

GREELEY 

1.23 

18.40 

78.  7 

CHARTREUSE 

0.67 

19.20 

94.3 

SCOTCH 

1.00 

2.04 

67.1 

KNICKERBOCKER 

0.63 

23.75 

96.7 

DURYEA 

0.55 

25.21 

90.8 

BOXCAR 

1.32 

10.60 

67.1 

REX 

0.67 

21.60 

70.0 

BENHAM 

1.40 

19.66 

121.0 

DUMONT 

0.67 

16.90 

79.9 

PILEDRIVER 

0.46 

12.06 

71.6 

PIRANHA 

0.56 

12.27 

77.2 

BRONZE 

0.53 

12.22 

50.7 

AGILE 

0.76 

10.97 

66.5 

WAGTAIL 

0.75 

14.59 

51.1 

FORE 

D .  49 

10.38 

23.0 

8UFF 

0.50 

16.37 

60.5 

NASH 

0.37 

16.26 

65.1 

TURF 

0.51 

20.65 

50.8 

CUP 

0.55 

11.13 

59.3 

KLICKITAT 

0.50 

20.40 

63.1 

BOURBON 

0.56 

11.60 

53.8 

AUK 

0.46 

30.47 

86.2 

PINSTRIPE 

0.30 

58.49 

1 14.  9 

PAR 

0.41 

5.80 

33.4 

COMMODORE 

0.76 

14.29 

69.0 

TAN 

0.56 

11.20 

63.8 

8ILBY 

0.  71 

10.76 

55.1 

CORDUROY 

0.69 

11.82 

63.3 

AVERAGE 

16.75 

69.4 

Res  tra i ned 

Shift 

Denth 

Error* 

Res  tra i ned 

a  (Sec) 

(km) 

(  km) 

o  (Sec) 

1.044 

7.05 

24.0 

0.4D4 

1 .094 

2.20 

7.3 

0.309 

1.056 

1  .  58 

17.5 

0.384 

1.244 

6.65 

35.3 

0.452 

1.065 

4.56 

26.1 

0.354 

0.997 

3.77 

22.1 

0.295 

0.731 

1.70 

17.9 

0.402 

1.312 

1.13 

2D. 7 

0.324 

1.017 

1.36 

4.4 

0.219 

1  .D25 

2.01 

0.4 

0. 33D 

1.048 

0.36 

8.2 

0.412 

O.L-16 

0.51 

(3.2) 

0.363 

1.147 

1 .  30 

8.2 

0.445 

0.552 

0.6D 

16.4 

0.298 

0.487 

5.85 

(24.7) 

0.428 

0. 791 

2.40 

7.9 

0.352 

0.948 

5.24 

4.7 

0.392 

0.652 

3.60 

11.4 

0.379 

0.700 

1.66 

20.5 

0.292 

0.661 

8.01 

39.9 

0.364 

0.714 

1.93 

(20.8) 

0.425 

0.693 

17.  36 

42.3 

0.352 

1.036 

11.33 

21.8 

D.252 

0-712  1.55  (29.0)  0.301 

1.009  r~  2.69  5.9  0.212- 

0-889  1.80  (6.7)  0.146 

0-959  1.23  0.4  0.137 

0.829  1.17  2.9  0.151 

0.902  3.90  14.9  0.355 


*  Depth  error  in  parentheses  indicates  that  solution 
zero  depths  were  used  in  computing  the  average. 


is  above  the  surface; 
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Concerning  the  depth  errors,  an  interesting  relation¬ 
ship  is  observed  between  the  standard  deviation,  o,  of 
least-squares  time  errors  of  the  depth  restrained  solution 
without  anomalies  and  the  depth  error,  dz,  of  the  corre¬ 
sponding  unrestrained  solution.  Figure  6  shows  the  results 
(dot-symbols)  and  the  linear  relationship  over  the  range  of 
standard  deviations  obtained,  which  can  be  expressed  as 


dz  =  z-zq  =  75a 

'  V 


where  zq  is  the  actual  depth  and  z  is  the  unrestrained 
depth.  The  relation  simply  states  that,  for  explosions  at 
NTS,  the  larger  the  least-scjuares  time  errors  (due  either  to 
anomalies,  to  reading  errors,  or  to  both),  the  greater  is 
the  depth  bias  in  the  unrestrained  solutions.  In  Figure  6, 
the  apparent  outliers  (denoted  with  arrows)  are  the  four 
events  with  the  least  number  of  stations:  11,  9,  11,  and  11 
in  order  of  increasing  a.  The  next  least  number  of  stations 
is  16  and  is  not  an  outlier. 

When  anomalies  are  applied  to  the  depth-free  solutions, 
the  results  are  as  shown  in  Figure  6  (x-symbols)  which  are 
clearly  below  the  linear  fit  to  the  no-anomaly  results. 
Therefore,  the  anomalies  apparently  remove  the  "depth-bias" 
(the  functional  relationshijxgbetween  a  and  dz)  for  NTS. 

The  remaining  15  km  average  "error,  corresponding  to  the 
3  kfn  average  error  in  epicenter,  is  caused  by  the  random 
errors  in  the  calibration  anomalies ,  and  by  the  random 
errors  in  arrival  times.  Also,  in  the  calculation  of  the 
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average,  all  negative  depths  were  set  to  zero,  although 
this  would  account  for  only  about  20%  of  the  remaining 
15  km  bias. 
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LOCATION  BIAS  EFFECTS 

In  the  preceding  analyses,  known  explosions  comprised 
the  data  set,  and  errors  caused  by  an  invalid  earth  model 
^re  eliminated  by  using  anomalies.  One  can  then  properly 
speak  of  location  accuracy  rather  than  precision  (relative 
accuracy).  In  the  case  of  earthquakes,  there  is  an  unknown 
error  m  the  assumed  position,  and  studies  concerning 
locations  using  earthquake  anomalies  are  always  suspect. 
However,  the  question  of  the  precision  achievable  by  applying 
anomalies  which  are  determined  from  an  event  whose  location 
is  unknown  can  be  addressed  if  a  calibration  explosion  is 
deliberately  mi 3 located. 

To  establish  a  base  from  which  a  comparison  can  be 
made,  true  anomalies  can  be  determined  from  the  actual 
position  of  BILBY  for  a  teleseismic  network  of  33  stations 
(Table  II)  and  applied  as  corrections  to  other  events 
occurring  within  the  same  anomaly  region.  Shifts  for  this 
case  calculated  with  depth  restrained  are  plotted  in  Figure  7 
for  the  27  events  (other  than  BILBY)  listed  in  Table  I. 

The  true  location  of  BILBY  is  indicated  on  the  figure. 

The  average  location  error  is  about  5  km  when  using  only 
the  BILBY  anomalies.  (The  normalized  plots  shown  in  Figure  7 
suggest  that  there  is  some  anomaly  instability  between  the 
two  sub-regions  and  that  one  should  probably  use  slightly 
different  anomalies  when  locating  events  in  subregion  II. 
However,  as  shown  in  a  previous  section,  if  an  anomaly 
average  of  several  events  is  used  instead  of  an  anomaly 
determined  from  a  single  event,  the  size  of  the  anomaly 
region  can  be  increased  without  sacrificing  accuracy. ) 
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If  the  actual  location  of  BILBY  is  now  deliberately 
given  an  error  of  +1°  longitude  and  +1°  latitude  (about 
140  km  to  the  northeast),  anomalies  can  be  re-determined 
from  this  assumed-correct  position  and  applied  as  corrections 
to  subsequent  events.  The  station  anomalies  that  one  would 
obtain  from  the  displaced  BILBY  location  are  given  in  Table  IX; 
teleseismic  anomalies  this  large  have  never  been  observed, 
so  it  is  safe  to  assume  that  bias  effects  will  generally  be 
smaller  than  those  imposed  on  this  example. 

The  results  of  locating  events  with  the  displaced-BILBY 
anomalies  are  shown  in  Figure  8.  The  average  location 
precision  (3.9  km)  is  about  the  same  as  the  average  location 
accuracy  in  Figure  7  (4.6  km).  (The  precision  is  actually 
a  little  better  than  the  accuracy  because  the  northeasterly 
bias  given  to  BILBV  artificially  moves  the  event  set  generally 
closer  to  the  networks,  thereby  increasing  the  sensitivity 
(dT/dA)  of  the  travel-time  tables.  Thus  the  solutions 
need  to  move  away  from  their  true  values  by  a  smaller  amount 
to  achieve  an  equivalent  minimization  of  the  residual  time 
errors).  The  main  point  to  be  made  is  that  regardless  of 
where  the  calibration  event  is  assumed  to  be  located,  the 
accuracy  of  subsequent  events  relative  to  the  calibration 
event  is  virtually  unaffected  if  anomalies  are  used.  This 
is  important,  for  example,  in  studies  of  nuclear  test  sites, 
of  event  migration  in  aftershock  sequences,  or  of  seismo- 
geologic  correlation  with  island  arcs,  oceanic  ridges,  and 
major  fault  zones  (Isacks  et  al  (1967);  Sykes  et  al  (1969); 
Mitronovas  et  al  (1969);  and  others),  in  which  the  pattern 
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TABLE 

IX 

j 

DISPLACED-BILBY 

ANOMALIES 

,  '  1 

Stati  on 
( C&GS ) 

.9 

Anomaly 
(Sec  )* 

Station 

(VELA) 

I 

| 

Anomaly 

(Sec)* 

ARE 

-13.52 

> 

.  BL-WV 

-4.77 

BEC 

-6.10 

BR-PA 

-4.10 

CAR 

-10.02 

CPSO 

-4.91 

FLO 

-1.40 

DH-NY 

-4.37 

6DH 

-3.62 

EB-MT 

1.11  ) 

GUA 

-18.10 

EU-AL 

-4.85 

KEV 

-7.90 

EU2AL 

-4.85** 

KIP 

-20.44 

GG-GR 

-6.40 

LPB 

-14.27 

HN-ME 

-3.44 

LPS 

-12.38 

LV-LA 

-4.99 

NNA 

-14.81 

NP-NT 

1 

-4.70 

NUR 

-7.54 

OO-NW 

-6.75 

DXF 

-4.10 

PZ-PR 

-8.21 

PTO 

-7.48 

RK-ON 

1 

0 

SCP 

-4.47 

SHA 

-5.23 

f 

1 

STU 

-7.04 

TOL 

-6.86 

TRN 

-12.57 

'  ; 

VAL 

-6.83  , 

I 

*  All  anomalies  are  referenced 

to 

station  RK-ON. 

**  Used  the 

anomaly  determined 

for 

EU-AL. 

;v- 


I, 

M' 
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of  events  (or  lack  of  it)  is  contributory  to  the  conclusions 
reached  regarding  the  location  of  clandestine  underground 
explosions,  seismicity,  or  global  tectonics. 
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LOCATIONS  USING  CONSTANT  NETWORKS 
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In  the  previous  sections  of  this  report,  the  locations 
were  made  using  every  available  recording  station,  the 
number  of  which  varied  between  9  and  49.  The  results 
obtained  without  anomalies  (Figure  2)  were  seen  to  be 
variable  in  magnitude  and  direction,  whereas  with  anomalies 
(Figure  4)  the  variability  was  significantly  reduced.  It 
will  now  be  shown  that  by  using  networks  composed  of  the 
same  stations  for  a  set  of  explosions  the  results  obtained 
without  anomalies  are  no  longer  variable  but  instead  yield 
exactly  the  same  precision  as  the  results  obtained  with 
anomalies.  Furthermore,  it  wi^l  be  shown  that  location  bias 
is  a  function  of  the  particular  network  used,  and,  because 
the  several  networks  investigated  have  essentially  identical 
quadrantal  distributions,  that  location  bias  cannot  be 
attributable  to  source  region  geology  alone  but  can  be  more 
properly  understood  as  a  function  of  the  total  travel 
path,  including  the  lower  mantle  and  the  crustal  and  upper  * 
mantle  geology  in  the  vicinity  of  the  stations. 

First,  a  five-station  network  is  selected:  RK-ON,  CPSO, 
COL,  HN-ME,  and  NP-NT.  This  network  has  an  azimuth  aperture 
of  109°  and  a  reasonable  range  in  distance  (21°  to  39°). 

Using  no  anomalies,  the  results  of  locating  17  events  recorded 
by  these  five  stations  are  shown  in  Figure  9  plotted  from  a 
common  origin.  The  average  error  is  22.9  km  and  all  of  the 
events  shift  northeast.  Clearly.,  this  particular  network 
may  be  said  to  have  a  northeast  location  bias  of  about  23  km. 

Using  the  same  network  but  with  the  NTS  anomalies  as 
given  in  Table  VI,  the  average  location  error  is  reduced  to 
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3.5  km  and  the  directions  are  randomized.  These  results 
are  plotted  in  Figure  10  to  the  same  scale  as  in  Figure  9.  The 
vectors  from  Figure  9  are  replotted  on  Figure  10  as  blue  lines; 
it  is  seen  that  the  relative  locations  (vector  terminal  points) 
are  identical  in  both  cases. 


If  to  the  five-station  network  used  above,  a  sixth 
station  (ARE)  is  added,  the  azimuth  aperture  is  increased 
to  157°.  The  results  without  anomalies  are  shown  in 
Figure  11a;  the  average  error  is  13.0  km  and  the  events 
again  shift  northeasterly.  Hence,  for  this  six-station 
network,  the  bias  magnitude  is  13  km,  and  the  bias 
direction  is  unchanged.  The  results  obtained  by  adding 
stations  one  at  a  time  and  without  anomalies  are  shown 
in  Figures  lib,  11c,  and  lid,  in  each  case  the  aperture 
remaining  at  157°.  The  average  location  errors  in 
Figures  9  and  11  can  be  summarized  as  follows:' 


$ 


Without  Anomalies 


Network 

Error,  km 

Direction 

No.  of 
Events 

Figure 

5-station 

22.9 

Northeast 

17 

9 

6-station 

13.0 

Northeast 

15 

•  11a 

• 

7-station 

10.8 

Northeast 

15 

lib 

8-station 

8.7 

Northeast 

,  12 

11c 

• 

9-station 

l 

9.8 

Northeast 

,  / 

11 

lid 
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Figure  11c.  Restrained  location  shifts  using  constant  eight- 
station  networks  and  no  anomalies  plotted  from 
common  origins. 
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STATIONS  USED  (9) 


AVCAAGC  CKKOR  •  9  §  Km 
H*S‘  9  0  Km  AT  if*  AZIMUTH 


*  «  « 
EAST  — - 


NAME 

distance 

*  azimuth* 

RK-ON 

21  0* 

430* 

OXF 

22  0* 

tt  4* 

CF30 

24  f 

•4  r 

COL 

S3  2* 

33§  0* 

NN-ME 

34  7* 

103* 

NP  -NT 

St  0* 

3St  f 

CAR 

31  •• 

IOt  0* 

KTO 

St  7* 

23  7* 

ARE 

it  2* 

133  2* 

RANGE  47. 2* 

.aperture 

•distance  ano  azimuth 

137  2* 

FROM  EVENT  1 

Figure  lid.  Restrained  location  shifts  using  constant  nine- 
station  networks  and  no  anomalies  plotted  from 
common  origins. 
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V  The  location  results  obtained  with  these  networks, 
but  with  anomalies  applied,  are  shown  in  Figure  12  with 
the  Figure  11  results  included  as  blue  lines,  again 
demonstrating  that  the  same  precision  is  achieved  if 
constant  networks  are  used.  The  results  of  Figures  10 
and  12  are: 

With  Anomalies 


Network 

Error ,  km 

Direction 

No.  of  - 
Events 

Figure 

5-station 

3.5 

Random 

17 

10 

6-station 

2.8 

Random 

15 

12a 

7-station 

2.8 

Random 

15 

12b 

8-station 

2.2 

Random 

12 

12c 

9-station 

1.9 

/  i 

Random  . 

11 

12d 

Therefore,  if  a  constant  network  is  used  to  locate 
events  from  a  particular  region,  anomalies  need  not  be 
applied  to  achieve  precise, locations.  However,  the  imposition 
of  a  constant-network  requirement  in  routine,  operation  is 
prohibitive;  stations  record  events  from  a  region  in  a  highly 
variable  way  due  to  their  detection  thresholds,  temporal 
microseismic  noise  fields,  local  seismic  activity,  and 
unpredictable  malfunctions  of  instruments.  It  is  much 
simpler  in  the  long  run  to  routinely  apply  anomalies  such 
that  the  locations  made  with  variable  networks  are  consistent. 

The  northeast- bias  in  Figures  9  and  11  is  not  to  be 
construed  as  a  function  of  the  NTS  regiop.  Alternatively, 
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Figure  12b.  Restrained  location  shifts  using  constant  seven 
station  networks  and  anomalies  plotted  from 
common  origins. 
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another  network  could  have  been  chosen,  using  stations 
FLO,  OXF,  COL,  HN-ME,  GDH,  and  TRN.  The  azimuth  aperture 
"this  network  is  127°,  and  the  stations  are  distributed 
azimuthally  from  north-northwest  through  east  to  southeast, 
as  they  were  for  the  networks  of  Figures  9  through  12. 

The  location  results  without  anomalies  are  shown  in 
Figure  13;  the  location  bias  for  this  network  is  nearly 
opposite  to  the  others:  51  km  to  the  southwest.  When 
anomalies  are  applied,  the  location  shifts  are  as  shown 
in  Figure  14  with  the  Figure  13  results  included  as  blue 
lines.  The  complete  results  for  this  six-station  network 


Without  Anomalies 

•6 

With  Anomalies 


Error,  km 
51'.  2 
5.57 


Direction 

Southwest 

Random 


No.  of 
Events 


Figure 


Similar  networks  could  be  selected  to  display  virtually 
any  bias  magnitude  and  direction  one  desires.  The  lack 
of  an  apparent  common  bias  among  similarly-distributed 
networks  definitely  eliminates  the  source  region  as  a  principal 
factor  in  causing  the  bias. 

Clearly,  bias  is  a  function  of  the  particular  network; 
more  precisely,  it  is  a  function  of  the  stations  uniquely 
defining  that  network.  But  bias  cannot  be  attributed  only 
to  "station  effect",  because  the  anomaly  would  then  be  a 
.  constant  (late  or  early)  or  nearly  so,  for  virtually  all 
regions.  Of  course,  some  azimuthal  dependence  might  be 
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observed  if  strong  dips  existed  under  the  station,  but 
a  large  proportion  of  the  anomaly  would  be  a  simple  DC 
shift.  It  is  well-established,  however,  that  the  anomalies 
are  not  constant,  are  not  simple  functions  of  azimuth,  and 
are  highly  variable  with  distance  as  well  (Chiburis  and 
Dean,  1965;  Chiburis,  '966b,  1968a).  The  variations  referred 
to  are  those  occurring  at  a  particular  station  for  many  event 
azimuths  and  distances  and  are  not  those  occurring  for  a 
particular  event  at  many  station  azimuths  and  distances 
(Cleary  and  Hales,  1966).  The  distinction  is  not  subtle. 

In  the  latter  case,  a  certain  amount  of  residual  averaging 
is  done  within  specified  windows  in  order  to  perceive  any 
unctional  variations.  Cleary  and  Hales  (1966)  have  found 
a  sinusoidal  variation  with  distance,  but  in  determining 
anomalies  for  location  purposes  (not  the  intent  of  Cleary 
and  Hales),  little  is  gained  unless  they  are  reproducible 
and  known  to  within  random  reading  error  (about  0.25  sec) 
for  each  region.  If  many  stations  are  contained  within 
a  window  of  either  distance  or  azimuth,  the  range  in  anomalies 
will  commonly  be  greater  than  2  sec  and  is  frequently  as 
high  as  5  sec,  certainly  too  large  for  predicting  station 
anomalies  for  accurate  location,  and  too  large  to  reasonably 
explain  all  bias  as  being  due  to  a  singl,  phenomenon. 

When  it  is  desired  to  determine  corrections  for  location 
purposes,  the  anomalies  at  a  particular  station  should  only  be 
averaged  within  an  anomaly  region,  because  from  one  region  to 
another  the  anomalies  change  substantially  (although  not 
erratically;  see  in  particular,  Chiburis,  1966a).  The 
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implication  is  that  the  major  portion  of  an  anomaly 
(especially  the  variation)  is  caused  by  slight  lateral 
and  vertical  inhomogeneities  within  the  mantle  between  the 
source  and  receiver,  the  effects  of  which  are  integrated 
along  the  entire  raypath.  Furthermore,  any  component  in 
the  anomaly  due  to  a  velocity  distribution  in  the  immediate 
vicinity  of  the  station  which  does  not  conform  to  the 
^-r’ave i— t ime  model  will  be  constant.  However,  the  observed 
distance  and  azimuth  variations  in  the  station  anomaly  are 
so  large  (5-6  sec,  in  many  instances)  that  they  probably 
mask  out  the  constant  term  entirely. 
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■  .  ANOMALIES  FROM  OTHER  EXPLOSION  REGIONS 

sio„/h<!  ^"°MlieS  giVe"  ^  TaMe  W  Mere  from  explo¬ 

it  a1"  a  1P0rtlOn  °f  NTS  Shom  in  Ti^e  1.  It  is  appropriate 
mined  f  ^T"  ^  thSSe  an°"alie=  and  those  deten¬ 
tion  f°"  eXF  °”S  ln  °ther  regtons .  Figure  IS  shows  the  posi- 

taons  of  seuerai  explosions  detonated  in  North  America.  Figure  16 

the  first  part  of  Table  X  give,  for  the  stations  used  in  this 
report,  a  comparison  between  the  NTS  average  anomalies  and  the 
anomalies  computed  from  the  other  explosions;  the  latter  part  of 
Table  X  includes  the  comparisons  for  an  additional  122  stations 
Whose  anomalies  were  determined  from  published  arrival  times. 

the  NTSth  thT  an°maly  COmpariso"s  “d  "ith  locations  made  using 
anomalies,  several  of  the  points  made  previously  in  this 

report  will  be  further  demonstrated  and  substantiated. 

Dlot/d  FlgUre  16’  ^  NTS  aVerage  anomaly  for  each  station  is 
ed  against  the  anomaly  determined  from  each  of  the  other 
explosions.  As  in  previous  tahl^c  t  • 

.  '  .  tables,  the  anomalies  are  with  respect 

NTS^ncTth  °f  C°UrSe’  anomalies  identical  for  both  the 

t  f  be  °ther  e>I,losio"s  "°uld  generate  straight  lines  on  this 
yp  of  plot  with  slopes  equal  to  one.  The  lines  shown,  however 
were  determined  from  a  least-squares  regression  of  the  La  in 
each  plot.  The  numerical  values  for  the' computed  slopes,  m,  are 
own  in  each  sub-figure.  The  standard  deviation  about  the  least 
squares  itted  line  is  also  presented  as  sy.x.  A  probable  caul 

line  not  going  through  the  point  CO,  0)  could  be  due  to 
an  anomalous  reading  at  RK-ON. 

The  event  JORUM  is  within  the  NTS  area;  hence,  as  expected 
e  s  ope  of  tb  fitted  line  for  the  JORUM/average  NTS  data, 

deviation  of ’ll  a  ^  ^  1’°  *  relatival>'  a”aH  standard 

n  of  fit.  At  the  other  extreme  are  the  results  for  the 
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Figure  16. 


Comparison  of  anomalies  between  the  Nevada  Test 
Site  and  other  regions. 


-53- 


*  Anoaaljr  For  PG2BC 

*  Anooaly  For  PG-8C 

••  All  Anomalies 

Referenced  To  Station 
RK-ON 

*  Anoaaly  For  S'F33B 

*  Anooaly  For  S*Z0B 


TABLE  X  (CONT'D.).  COMPARISON  OF  ANOMALIES  FOR  SELECTED  EXPLOSIONS 


V 


% 


event  LONG  SHOT,  located  in  the  Aleutian  Islands  with  a  slope  of 
0.52  and  a  standard  deviation  of  slightly  greater  than  1  sec, 
(Figure  16H).  For  the  event  FAULTLESS,  which  is  approximately 
100  km  north  of  the  NTS,  the  slope  is  0.88  with  .a  standard  devia¬ 
tion  0.69  sec,  which  indicates  a  greater  disagreement  with  the 
NTS  anomalies  than  was  evidenced  in  the  JORUM  anomalies.  Each  of 
the  other  explosions  also  show  definite  variations  from  the  NTS 
values ,  either  in  the  slope  of  the  fitted  line  or  in  the  standard 
deviation  of  the  fit,  or  both.  The  standard  deviation  of  fit  for 
the  GNOME  explosion  is  only  0.26  sec  but  the  slope  is  0.72,  which 
is .far  enough  from  1.0  to  cause  a  large  error  when  the  average 
NTS  anomalies  are  used  for  locating  GNOME.  The  comparison  for  the 
RULISON  anomalies,  with  a  slope  of  .96  indicates  good  agreement 
on  the  average  with  the  NTS  anomalies,  but  the  standard  deviation 

-*-s  1*^3  sec,  which  is  too  large  to  expect  good  location 
results.  '  ''  • 

For  many  of  these  events,  the  epicenter- to-station  azimuths 
are  similar,  but  the  anomalies  are  significantly  different;  the 
differences  are  so  great  that  the  cause  of  the  anomalies  cannot 
reasbnably  be  attributed  to  local  geology  ht  the  recording  sta¬ 
tion.  For  example,  NP-NT  has  an  anomaly  range  of  1.73  seconds  and 
KEV  3.53  seconds,  changes  which  can  hardly  be  due  to  station  effect, 

The  largest  anomaly  variation  for  any  one  of  tiie  stations 
listed  in  the  first  part  of  Table  X  is  -5.55  seconds  at  AAM;  for 
any  one  of  the  stations  in  the  entire  table,  6.28  seconds  at  BCN. 

The  largest  range  in  the  anomalies  for  a  particular  event  is  6.65 
seconds  (-0.37  sec  to  +6.28  sec)  for  SALMON.  The  largest  range  for 
any  station  and  any  event  is  9.43  seconds  -C^3. 15  sec  to  +6.28  sec). 

Locations  of  the  events  with  deaths  restWined-and  using 
every  available  station  in  Table  X,  but  withouV" anomalies ,  are 
summarized  in  Table  XI.  If  it  is  assumed  that  the  NTS  anomalies 
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are  valid,  locations  made  using  them  can  be  compared  to  the 
no-anomaly  solution  tc  note  the  improvement  in  location  accuracy. 

a  les  XII  and  XIII  respectively  give  the  no-anomaly  and  anomaly 
lecatton  results  obtained  using  only  those  stations  with  a  listed 
NTS  average  anomaly.  Comparing  the  two  tables,  essentially  no 
improvement  is  achieved  by  using  NTS  anomalies;  for  some  events 

rnp|1rSe  diStanC8S  fr°m  NTS>  locations-  are. worse.  Excluding 
,  w  ich  is  in  the  NTS  area,  the  average  location  error  for 
the  seven  events  obtained  without  anomalies  is  -20.2  km  and  with 
anomalies  is  20.8  km.  The  ratio  of  the  location  error  obtained 
without  th  NTS  anomalies  to  that  obtained  with  the' NTS,  anomalies 
is  presented  in  Table  XIV.  As  expected,  JORUM  shows  by  far  the 
greater  improvement  in  location  accuracy.  Also  presented  are  the 
ratios  of  the  variances  of  time  errors  from  the  least-squares 
ocations  for  each  explosion.  Again,  the  ratio  for  JORUM  is  sig¬ 
nificant,  and  so  too  is  the  ratio  for  FAULTLESS.  The  ratio  for 
GNOME,  7.4,  although  very  large  is  not  statistically  significant 
using  the  "F"  test  because  there  were  only  five  stations  used  in 
the  location  scheme.  The  only  other  explosion  with  a  reasonable 
improvement  in  variance  was  SALMON  but  again,  too "few  stations 
were  used  for  the  ratio  to  be  significant. 

Therefore,  it  is  known  that  the  anomalies  change  signifi- 
cantiy  for  most  stations  from  one  region  to  another  nearby,  and 
at  for  all  stations  they  can  change  when  the  regions  are  far 
apart.  Consequently,  if  for  no  other  reason  the-  that  the  labor 
expended  in  computing  useless  anomalies  is  saved  and  that  the 
location  results  are  less  apt  to  be  misinterpreted,  it  is  better- 

location  procedures  to  apply  no  anomalies  at  all  than 
the  wrong  ones. 
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CONCLUSIONS 


A  location  study  was  made  of  28  underground  explosions 
detonated  in  the  northern  area  of  the  Nevada  Test  Site 
(NTS).  The  events  were  large  enough  such  that  ?~wave 
tinning  errors  were  within  acceptable  limits. 

The  networks  used  for  locating  the  events  were  comprised 
of  between  9  and  49  teleseismic  stations  having  two-  to 
three-quadrant  distributions.  Because  all  of  the  events 
have  known  positions,  the  results  regarding  the  effectiveness 
of  travel- time  anomalies  on  location  accuracy,  location 
precision,  depth  error,  etc.,  are  more  consistent  and  simpler 
to  evaluate  than  if  the  events  were  earthquakes  of  unknown 
positions . 

The  conclusions  inferred  from  this  study  are  as  follows: 

1.  Depth-restrained  location  errors,  using  all 
available  stations  without  anomaly  corrections, 
average  about  7  km  but  are  as  large  as  20  km. 

The  lack  of  consistency  among  the  locations 
suggests  that  the  errors  are  not  due  solely  to 
conditions  in  the  vicinity  of  the  source  region 
(NTS).  Also,  if  different  travel-time  tables 
are  used,  the  average  error  is  unchanged,  but 

a  significantly  different  pattern  of  locations 
is  obtained. 

2.  Depth-restrained  location  errors,  again  using 
all  available  stations  but  with  anomaly 
corrections  determined  from  a  few  calibration 
events,  are  reduced  to  2. 5-3.0  km.  Comparing 
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these  results  with  those  of  a  previous  study 
(Chiburis,  1968),  the  same  average  error  of 
2. 5-3.0  km  is  achievable  regardless  of  the 
network  used.  In  the  present  study,  two- 
and  three-quadrant  networks  are  used,  and  the 
average  error  without  anomalies  is  7  km.  In 
the  previous  study ,  many  of  the  networks  were 
single  quadrant,  and  the  average  error  with¬ 
out  anomalies  was  about  25  km.  However,  when 
anomalies  are  applied  the  average  error 
remains  the  same,  and  network  stability  is 
less  a  factor  in  producing  location  errors. 

The  size  of  the  area  at  NTS  across  which  the 
anomalies  are  valid  is  at  least  25  km  by  70  km 
as  evidenced  by  the  comparable  location  errors 
obtained  at  separated  sub-regions  within  NTS. 
This  implies  that  the  task  of  calibrating 
other  specific  regions  of  interest  would  not  be 
formidable . 

If  average  anomalies  are  used,  instead  of 
those  from  a  few  calibration  events,  the 
location  errors  can  be  reduced  to  approxi¬ 
mately  2  km.  This  represents  the  best 
possible  teleseismic  accuracy  for  NTS  events. 

Unrestrained-depth  solutions,  without  applying 
anomalies,  yield  depth  errors  of  about  70  km 
and  epicenter  errors  of  about  17  km.  All  of 


the  depth  errors  are  positive  (solutions 
too  deep)  for  the  networks  involved. 

A  linear  relationship  is  observed  between  the  least- 
squares  standard  deviation  a  of  time  errors  of  a  restrained 
solution  and  the  depth  error  dz  of  the  corresponding 
unrestrained  solution:  dz(km)  =  75 a  (sec). 

/ 

6.  Unrestrained-depth  solutions  obtained  with 
the  application  of  anomalies  yield  depth 
errors  of  about  15  km  and  epicenter  errors 
of  4  km,  improvement  factors  of  about 
five  and  four,  respectively,  over  the  no¬ 
anomaly  solution. 

A  linear  relationship  is  no  longer  observed  between  a 
and  dz,  the  relation  now  appearing  random;  hence,  the 
anomalies  effectively  remove  the  NTS  "depth  bias". 

7.  Relative  accuracy  (precision)  is  unchanged, 
even  if  the  calibration  events  are 
deliberately  mislocated  by  140  km;  that  is, 
if  anomalies  are  determined  from  an  event 
having  an  assumed-correct  location,  the 
location  errors  of  subsequent  events 
relative  to  the  calibration  event  -remains 
2. 5-3.0  km  with  the  anomalies  applied, 
regardless  where  the  calibration  event 
actually  is. 

8.  When  using  a  network  composed  of  the  same 
stations  for  an  entire  set  of  events,  the 


-66- 


r 


*'**V'-im*.n***!r, 


W 

H  • 

3* 


resultant  "bias"  (magnitude  and  direction 
of  the  location  shifts)  is  a  consistent 
and  unique  function  of  that  network;  when 
i  another  constant  network  is  used  to  locate 

the  same  events,  an  entirely  different  "bias" 
emerges.  The  apparent  lack  of  a  common  bias 
eliminates  the  source  region  as  the  principal 
cause  of  the  anomalies . 

Location  bias  may  then  be  considered  in  large  measure 
as  the  result  of  slight  lateral  and  vertical  inhomogeneities 
within  the  mantle  between  the  source  and  receiver,  the 
effects  of  which  are  integrated  along  the  entire  path.'"' 

9.  If  a  constant  network  is  used,  the  relative 
locations  obtained  without  anomalies  are 
identical  to  those  obtained  with  anomalies, 
except  for  a  bias  translation  appropriate 
for  that  network.  Therefore,  anomalies 
need  not  be  applied  for  precise  location 
work,  if  the  requirement  of  always  using 
the  same  network  is  not  too  stringent. 

10.  Station  anomalies  determined  from  explosions 
occurring  in  regions  other  than  NTS  show 
large  differences  when  compared  to  NTS,  some 
.  exceeding  6  sec. 

When  using  the  NTS  anomalies  to  locate  these  other 
events  essentially  no  improvement  is  made  in  the  solutions: 
the  average  error  without  anomalies  for  seven  events  not 
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in  the  NTS  drea  is  20.2  kin  and  with  anomalies 

20.8  km.  In  general,  it  is  better  to  apply  no  anomalies. 

at  all  rather  than  the  wrong  ones. 
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